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Abstract 

We present techniques for identifying and analyzing galaxy groups and show that they provide a powerful 
and robust discriminator between cosmological models. We apply these methods to high-resolution particle- 
mesh (PM) N-body simulations of structure formation in three SI = 1 cosmological models; Cold plus Hot 
Dark Matter (CHDM) with £l co id = 0.6, £1„ = 0.3, and Slbaryon = 0.1 at b=1.5 (COBE normalization) 
and two CDM models; b=1.5 and b=1.0 (COBE normalization). Groups are identified with the adaptive 
friends-of-friends algorithm of Nolthenius (1993). Our most important conclusions are: (1) Properties of 
groups are a powerful and robust discriminator between these Gaussian cosmological models whose spatial 
and velocity dispersion properties differ on ~Mpc scales. We test robustness against several methods for 
assigning luminosity to dark matter halos, for merging CfAl data, and for breaking up massive dark matter 
halos to correct for the overmerger problem. (2) When allowance is made for the higher than typical large 
scale power present in the CfAl data, CHDM at our £1„ = 0.30 produces slightly too many groups and too 
high a fraction of galaxies in groups, while the fraction grouped in CDM is far too low. A slightly lower £1„ 

CPS would appear to produce excellent agreement with all measures, save one: for all simulations, median group 

sizes are up to a factor of 1.7 larger than equivalently selected CfAl groups. This is either a real difference 

Oh or perhaps due to residual resolution limitations. (3) The standard group M/L method applied to our SI = 1 

simulations gives SI ~ 0.1 for CHDM and CfAl, and SI ~ 0.35 for CDM. We show quantitatively how three 

-4_^ different effects conspire to produce this large discrepancy, and conclude that low observed Si's need not 

argue for a low SI universe. When overmergers are broken up, the median virial-to-DM mass M„; r / 'Mom 
of 3D selected groups is ~ 1 for all simulations. Groups with Mom > 10 14 Mq appear virialized in all 
simulations. We measure global velocity biases b v similar to previous studies. Within 3D-selected groups, 
CHDM and CDM b=1.5 show a stronger bias of b v = 0.7 - 0.8, while CDM b=1.0 shows group 6„'s - 1- 

Subject headings: cosmology: theory — dark matter — large scale structure of the universe — galaxies: 
formation — galaxies: clustering 



1. Introduction 

The simplest viable scenario of cosmic evolution begins with an SI = 1 expansion seeded from inflation 
with Gaussian primordial fluctuations with a Harrison-Zel'dovich scale invariant spectrum. These fluctua- 
tions later collapse via gravitational instability to form structure from galaxies on up. This picture has been 
remarkably resilient in the face of steadily mounting observations (Gorski, et al. 1994, Dekel 1994). The 
nature of the dark matter will govern how structure forms, and is thus in principle recoverable from the 
statistics of such structure (and hopefully, of course, from direct detection as well). Evidence of structure on 
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very large scales in the late '80 's began a series of efforts which have all but ruled out standard Cold Dark 
Matter (CDM, e.g. Baugh & Efstathiou 1993). This has motivated the investigation of alternative models 
with the desirable properties of significant power on large scales and small power on small scales, including 
Cold + Hot Dark Matter (CHDM). As recently emphasized by Pogosyan (1994), standard CDM's demise 
leaves CHDM as the only remaining theory compatible with the simplest and most aesthetic versions of 
inflation. Even more interesting, recent neutrino oscillation experiments provide preliminary evidence that 
at least one neutrino species indeed has a cosmologically important mass (Caldwell 1994, Primack et al. 
1994). 

Beginning with Klypin, et al. (1993, hereafter KHPR), we explore in this series of papers the conse- 
quences of a universe dominated by Cold + Hot Dark Matter and compare it with standard Cold Dark Matter 
and with observations. The CHDM model with £1„ = 0.3 has already shown good agreement with observa- 
tions of the galaxy correlation function (Baugh & Efstathiou 1993), galaxy pairwise velocities (Somerville et 
al. 1995) and bulk velocities (KHPR), the cluster-cluster correlation function (Holtzman & Primack 1993; 
Klypin & Rhee 1994), the variance and skewness of the Abell/ACO cluster distribution (Plionis, et al. 1994), 
the amplitude of the power spectrum from POTENT reconstruction of the local density field (Seljak & 
Bertschinger 1994), the QDOT-IRAS Redshift Survey power spectrum (Feldman, Kaiser, & Peacock 1994, 
Fisher et al. 1993), quietness of the local Hubble flow (Schlegel, et al. 1994), the x-ray properties of clusters 
vs. redshift (Bryan, et al. 1994) and an initial analysis of the properties of galaxy groups (Nolthenius, Klypin 
& Primack 1994 a; hereafter NKP94). More controversial, (Ghigna, et al. 1994) find the standard CDM 
void probability function is in good agreement with observations in the Perseus-Pisces Supercluster region 
while that for CHDM is too high. deLapparent, Geller & Huchra (1986), however, find standard CDM 
produces voids too small, relative to those in the CfAl Survey (Davis, et al. 1982), suggesting CHDM would 
be in better agreement here. Also, while CHDM agrees better with the APM galaxy angular correlation 
function oj(6) than does standard or tilted CDM, CHDM's u>(9 > 3°) is still slightly too low (Yepes, et al. 
1994). Finally, the Hubble constant Ho can be no larger than ~ 50 km s _1 Mpc _1 to avoid problems with 
overproducing clusters and, as for all SI = 1 models, staying within cosmic age constraints. Therefore, if 
the current observational evidence for Ho — 80 km s _1 Mpc _1 is confirmed, then CHDM in its simplest and 
most aesthetic form is ruled out. Throughout this paper, we assume Ho = 50 km s _1 Mpc _1 . Aside from 
this, the easiest way to rule out CHDM is to find massive collapsed objects at high redshift, as CHDM forms 
structure significantly later than most competing models. These observations are currently underway, and 
indeed the first results are pushing the preferred £1„ down to perhaps 0.2 (Klypin, et al. 1994a). 

In this paper, we show how the statistics of galaxy groups provide a powerful discriminator between 
cosmological models. We identify and analyze groups from the 100 Mpc 512 3 particle-mesh simulations 
described in Klypin, Nolthenius, & Primack 1994b (hereafter KNP). The box size was chosen to give a good 
statistical sample of galaxy groups. With a mesh size of 195 Kpc (perhaps typical of an ~ L* galaxy dark 
matter halo), we do not have the resolution to study individual galaxies. However, typical galaxy groups 
have virial radii of ~ 1 — 2 Mpc and begin to be resolved. In NKP94 we presented our methods and an 
initial analysis of groups in the CHDM and CDM simulations, and in the CfAl Survey. This paper provides 
a more complete analysis, and includes corrections for effects not considered earlier. 

It can be validly argued that the properties of real groups depend sensitively on the properties of 
real galaxies, since galaxies are still the tracers for which we have the best observational statistics, and 
so the problem of realistically identifying galaxies is not so easy to circumvent. For example, group virial 
radii depend strongly on the spatial distribution of individual galaxies. Until more realistic large-scale 
simulations become feasible, the only way to deal with this problem is to add in by hand properties which 
are well determined physically but poorly understood or difficult to simulate (e.g. luminosity), and to do so 



2 



equitably across all competing cosmological models. We have addressed the resolution problem by filtering 
the observational data to the same resolution as the simulations. Fortunately, tests on most of our chosen 
statistics show little sensitivity to the relatively poor spatial resolution and related uncertainties, as we show 
in the present paper. A more significant problem is that of "overmerging" (Katz & White 1993, Gelb & 
Bertschinger 1993). Dark matter (DM) halos are extended, soft objects which merge easily. Galaxies undergo 
dissipative collapse to denser, smaller-cross section objects which merge at perhaps only half the DM merger 
rate (Evrard, Summers, & Davis 1993, hereafter ESD). We considered several schemes for breaking up our 
most massive halos. While with reasonable assumptions group rms velocities and fractions grouped appear 
insensitive to the details of halo breakup, we will show it is nevertheless possible to find breakup prescriptions 
that do lead to significantly different group properties. 

We analyze groups selected from two different versions of our catalogs. In order to gain insight into 
the properties of our models using all available information, we select groups from the complete 100 Mpc 
simulation boxes using full three dimensional (3D) information ("box groups"). To make meaningful, direct 
comparisons with the CfAl data, we first make magnitude limited "sky catalogs" in redshift space from each 
simulation, and then identify groups using identical criteria. 

The construction of the 512 3 grid of cells within the 100 Mpc box and the simulation calculation methods 
are described in KNP. In NKP94, we presented a robust and discriminatory statistic, rms group velocity 
i> gr (the rms velocity of all galaxies within the group, i.e. the conventional "velocity disperion") vs. the 
fraction of catalog galaxies in groups / gr , and compared CHDM with CDM and with the CfAl Survey 
data. Here, we present a more complete description of the construction of improved galaxy redshift catalogs, 
including corrections for some effects not considered earlier, analyze the properties of 3D-selected groups 
from the full box, and present other comparisons between these catalogs and observations. An accompanying 
computer visualization video sequence (Brodbeck, et al. 1994) compares visually the differences between the 
simulations, sky catalogs, and the real universe. CDM structures are rather puffy filaments, with clusters at 
the intersections. By contrast, CHDM is seen to have remarkably delicate filaments which are reminiscent, in 
a picturesqe way, of a well prepared egg flower soup. Blurring in redshift space reduces the visual impression 
of these differences, yet their effect remains powerful on the statistics studied here. 

Grouping in redshift space is done using the adaptive algorithm of Nolthenius (1993; N93), although tests 
with the original algorithm of Nolthenius and White (1987; NW) gave essentially the same results but with 
slightly lower group velocity dispersions. Recently, Frederic (1995a, b) has claimed that the NW grouping 
algorithm seriously underestimates group velocity dispersions and hence M/L ratios when compared with 
simulation groups. However, these studies fail to emphasize the importance of the underlying cosmology 
to setting the grouping algorithm's appropriate link in redshift. Using only a low bias standard CDM 
simulation, as Frederic did, will indeed show that the NW redshift link normalization of V5 = 350 km s _1 
is too low, as already shown in NKP94, due to the high pairwise velocities in this model. The fault lies 
not in the algorithm, but in the assumption that standard CDM is an appropriate calibrating cosmology for 
determining the merits of grouping algorithms used on real observations. Likewise, discrepancies he finds 
with the Ramella, et al. (1989; RGH) CfA2 group results vanish when one realizes the CfA2 dataset should 
not be assumed to have CDM-like properties. 

2. Construction of the Galaxy Halo Catalogs 

A detailed description of the simulation calculations is given in NKP94 and KNP. Briefly, we use a 
particle mesh code on a 512 3 grid, with 256 3 cold and 2 x 256 3 hot particles. The cold particle mass is 
2.9 x 10 9 M Q and 4.1 x 10 9 M Q for CHDM and CDM, respectively. As before, we refer to the CDM b=l 
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(COBE normalization, see Smoot, et al. 1992) and b=1.5 simulations as CDM1 and CDM1.5, respectively. 
The two CHDM simulations, CHDM1 and CHDM2, are both at b=1.5 (COBE normalization), and differ 
only in their initial conditions. The CHDMi, CDM1, and CDM1.5 simulations all began with the same 
random number set describing the amplitudes of the initial waves perturbing the particles. It was later 
found that this random number set had abnormally high power on large scales. The power spectrum was a 
factor of ~ 2 higher than typical on scales comparable to the box size. The probability of this occuring in 
any given realization was estimated at ~ 10%. A second CHDM simulation, CHDM2, began with a much 
more typical spectrum. This fact will be important later in interpreting the comparisons with CfAl data. As 
it turns out, there is good evidence that CfAl has unusual large scale power as well. The value of beginning 
each competing model with the same random perturbation set is that it guarantees the same large scale 
structures will emerge in each, so that differences between simulations will solely be due to differences in the 
underlying physics of the evolution and not to cosmic variance. 

A galaxy halo (boldly referred to hereafter as a "galaxy") is defined as a mesh cell with a sufficiently 
high dark particle mass overdensity (6p/p) cu t at the end of the simulation (z=0). We quantify the mass of a 
halo by the mass contained within a cell (or sometimes 3x3x3 cell) boundary. This is obviously crude. The 
boundaries are arbitrary and no attempt is made (nor is it possible) to include only the gravitationally bound 
particles; our spatial and force resolution is too poor to justify such refinements. Nevertheless, as long as 
galaxy luminosity monotonically rises with 1-cell mass, our results are insensitive to how mass is quantified. 
When our purpose is to optimally delineate structure we keep all cells above (6p/p) cu t = 30, giving N=29151 
halos in CHDMI, N=37,164 in CDM1, N=45,592 in CDM1.5, and N=29,795 in CHDM 2 . When constructing 
magnitude limited redshift catalogs and attempting to match CfAl galaxy number densities, experiments 
showed that (6p/p) cu t = 80 — 150 was best, cutting galaxy totals by about 70%. Note that this is slightly 
less than the (6p/p) cu t = 170 corresponding to virialization, (e.g. Kaiser 1986), and seems reasonable if one 
assumes virialization should actually apply to a denser core of material closer to the visible galaxy. Since 
each cell is 1/512 of the 100 Mpc box, or 195 Kpc , a cell would most properly correspond to a dark matter 
halo surrounding a typical L* galaxy. Two measures of the mass of such cells were calculated; the dark 
particle mass contained within the cell, and the dark particle mass contained within a 3 3 cell cube (3-cell) 
centered on the cell of interest. 

3. Breaking Up Massive Halos 

One uncertainty in our earlier results is that the overmerging of dark matter galaxy halos may have 
significantly lowered the rms velocities of galaxies within groups, and perhaps also lowered the total fraction 
of galaxies in groups. The dominant galaxy will tend to sit near the group center and have low center- 
of-mass velocity. By contrast, baryonic dissipation causes earlier collapse into several smaller but higher 
velocity galaxies (Katz & White 1993) whose merging rate is slower (ESD). If, as has been argued by some, 
the overmerging seen in numerical simulations is mostly due to poor mass resolution, our results should be 
relatively secure. Our mass resolution is quite good; an L* galaxy has about 100 cold particles and twice as 
many hot particles. Assuming the more likely probability that the lack of dissipation is indeed the dominant 
source of overmerging, we now ask above what halo mass M bu does it become important? Katz & White 
(1993) find that some individual halos could be as massive as 1.5 x 10 14 M Q . However, using reasonable M/L 
ratios, cooling efficiencies, and star formation rates, they argue most halos above M bu = 5 x 10 13 M Q should 
be broken up, and that M bu could be as low as 2.8 x 10 12 Mq . Gelb (1992) finds too many halos above 
Vdrc = (^t^-) 1 / 2 = 350 km s _1 which, for r corresponding to our cell size, is equivalent to 3.4 x 10 12 Mq . 
ESD also address this issue, and with better mass resolution they find a lower M bu limit of 7 x 10 11 Mq . Our 
most massive CHDM halo has 3-cell mass of 5 x 10 13 Mq (or a 1-cell mass of 5 x 10 12 Mq ). To date, ESD is 
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the best published work available on galaxy formation in groups and clusters in a dark matter background. 
Their Figure 12a provides a rough relation between dark matter halo mass and the number of galaxy-like 
objects ("globs" in their nomenclature) within the halo in standard CDM. While this figure applies to z=l, 
their most recent results show the number of globs stays fairly constant to z=0, at least in standard CDM 
(Evrard 1994). They define the halo's mass as the mass within the sphere containing an average overdensity 
corresponding to virialization: 6p/p = 170. Assuming unevolving halo size from z=l to the present, the 
cosmological expansion factor of 2 then gives a corresponding hp/ p = 2 3 x 170 = 1360 for Mj, u in our z=0 
simulations. We found the radius r e jj and mass M e ff of this sphere by determining the radii of spheres 
encompassing the volume of our 1-cell and 3-cell masses, and assuming density fell linearly from the 1-cell 
sphere radius to the 3-cell sphere radius. The maximum radius allowed to enclose a single breakup candidate 
correponded to the sphere enclosing a 2 3 -cell volume, or halfway to the nearest allowable halo. This was done 
to avoid double counting some of the exterior mass. We then broke up halos whose mass M e ff within a sphere 
of average overdensity 6p/p = 1360 was greater than 7 x 10 11 Mq . The result of this procedure was N om = 
138 (CHDMi), 157 (CHDM 2 ), 571 (CDM1.5), and 658 (CDM1) halos identified as overmergers and suitable 
for breaking up into fragments. If mass is roughly proportional to light, observed galaxies suggest that 
fragments should be assigned Schechter-distributed masses. Dissipational hydrodynamic codes also produce 
gaseous galaxy-like objects with Schechter distributed masses (Evrard 1994). We therefore constrained the 
fragments to follow a Schechter distribution with the characteristic mass M* = (M/L)L* = 4.3 x 10 11 Mq 
(CHDM) or 1.08 x 10 12 M e (CDM). L* is 4.3 x 10 10 L e from the CfAl and (M/L) is the median M/L 
of the simulation halos; 10 for CHDM and 25 for CDM. The faint end slope was set to the merged CfAl's 
a = —1.26 (see §6). Each overmerger was replaced with fragments of total mass M, where M was the 
virialized overdensity mass M e jj calculated above, minus the mass expected to be in fragments below the 
overdensity limit (6p/p) cu t selected for the simulation (and thus too faint to see); 80 for both CHDM and 
156 for both CDM simulations (see §6). For all simulations, 16 — 19% of the integrated Schechter function 
is contained in fragments below (6p/p) cu t. For each overmerger, the brightest fragment was given a mass 
Mtj = (Mb u M e j j) 1 / 2 (Evrard 1994) and placed at the original overmerger's position. Remaining masses 
were then randomly selected from a Schechter distribution and randomly placed into any halo which could 
accept it without overfilling. This continued until the next random fragment mass was too large to be added 
to any halo. At this point, all halos were ~ 99% full of fragments. Each set of fragments for a given halo 
was then ordered by mass and then sequentially placed as close as possible to the parent halo cell while still 
enforcing the 2-cell closest neighbor resolution limit. Thus the most massive fragments were placed closest 
to the original DM center. Each fragment was then given a randomly oriented, random Gaussian velocity 
with dispersion equal to the rms velocity V°™ gh of all halos within 1 Mpc of the overmerged halo (or, in the 
few cases this did not include at least 4 halos, out to the 4-th nearest halo). The median values (V°™ h ) me d 
of these neighborhood rms velocities were ~ 300 km s _1 for all models, but individual overmergers ranged 
as high as ~ 1200 km s _1 (CHDM) to ~ 2000 km s _1 (CDM1). The most massive fragment was left at the 
original overmerger's position. Observations indicate that the brightest (assumed the most massive) galaxies 
in groups and poor clusters are moving at a center-of-mass velocity only ~ 0.25 that of their lower luminosity 
neighbors (Bird 1994). We therefore multiplied the brightest fragment's velocity by 0.25. The fragments 
of a given overmerger do not, at this point, satisfy conservation of momentum. In the frame of the parent 
overmerger, the net momentum of the N fragments with masses rrii is m-iv} = p. To enforce momentum 
conservation we correct each fragment velocity by adding a velocity differential 6v = —p(%2f m 8 ') _1 . Finally, 
the masses were rescaled by Mi ce u/M e ff so that they followed the 1-cell convention used for all other 
cataloged halos. 3-cell masses were found in a similar way. 
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The simulation halos we analyze below have all been subjected to this, our preferred breakup scheme. 
We also experimented with other breakup schemes. These are described in detail below, and summarized in 
Table 1. 

The first alternative scheme, hereafter "Method 1" made use of ESD's Figure 12a, giving the number of 
galaxies per DM halo vs. DM halo mass. Mj, u was again 7 x 10 11 Mq , but this time the halo mass was simply 
assumed to be the 1-cell mass; generally a bit lower than the mass inside a sphere of average overdensity 
1360. We broke each overmerger candidate into ESD's nominal number of fragments, giving them essentially 
equal masses (when ordered by mass, each fragment was an arbitrary 10% more massive than the previous 
fragment). Velocity assignments were as before, except rather than using all neighbors inside a distance of 
1 Mpc, we used the 10 nearest neighbors, whatever their distance (typically out to a distance of ~ 2 Mpc ; 
as large as a medium sized group). The most massive fragment's velocity was not reduced. Placements 
again enforced the 2-cell nearest neighbor limit, but this time without regard to putting the most massive 
fragment closest to the DM center. The largest overmerger in any simulation spawned only 7 fragments, and 
the large majority of overmergers produced only 2 or occasionally 3 fragments. By essentially maximizing 
the mass of each fragment, this prescription guaranteed (albeit unintentionally) the highest possible number 
of fragments surviving the magnitude limit and making it into the final sky catalogs (see §6). Since these 
galaxies survived in close pairs or groups, it also significantly raised the fraction of galaxies in groups. 

The second scheme ("Method 2") made the more reasonable assumption of Schechter-distributed masses, 
randomly selected and assigned with the same M* parameters as our adopted scheme, but with ESD's "glob" 
faint end slope of a = —1.35. Mi, u was the same as for our adopted scheme. We did not constrain the brightest 
fragment's mass or velocity (except by enforcing momentum conservation, as before). Fragments were 
positioned randomly, regardless of mass, but enforcing the 2-cell limit. By not insuring at least one reasonably 
massive fragment, the result was a larger number (~ 30 on average) of low mass fragments, relatively fewer of 
which ultimately survived the magnitude limit for inclusion into the sky catalogs. Velocities were randomly 
sampled from a Gaussian distribution with dispersion equal to the circular velocity V° m = (GM e j j /r e j j) 1 / 2 
associated with the virialization mass and radius described in our final method above. 

Our preferred scheme assigns fragments their nearest neighbors' rms velocity V°™ gh rather than the 
circular velocity just described, for several reasons. First, the most appropriate mass and radius to use 
are poorly defined. Second, galaxies clearly form early, and, in the dense environments common to these 
overmergers, it is the larger tidal field of the group which seems likely to eventually determine final fragment 
velocities. Finally, one of our goals was to evaluate how well the virial theorem measures mass in these 
systems, and enforcing dynamically determined velocities would bias these results. Figure 1 shows a scatter 
plot of V c om vs. V°™ gh for each overmerger in the CHDM and CDM simulation boxes. The V c om and V°™ gh 
distributions differ substantially. However, their median values are very similar (V° m ) me d — (V,?™gh)med, as 
shown in Table 2. Circular velocities are slightly higher than V°™ gh ; by ~ 5% for CHDM and by ~ 15% 
for CDM. Thus, using circular velocities would likely have raised our final group rms velocities by only a 
few percent, more so for CDM than for CHDM. (However, it might also be argued that a more appropriate 
dynamical velocity would be lower than V c , closer to the virial velocity, which would be lower by ~ 2 -1 / 2 ). 
Note also that the extended tail of high V°™ gh is especially pronounced for the CDM models, and is related 
to CDM having higher pairwise velocities than CHDM (KHPR). 

Another approach ("Method 3") is to leave the simulation halos alone and instead attempt to put real 
galaxies back into "overmerged" halos by merging them, i.e. taking a luminosity weighted averaged of their 
positions and redshifts, then combining luminosities. In this case, one wants to merge CfAl galaxies which 
are within about 1.5 — 2 cells (~ 350 Kpc ) on the sky, and within a velocity separation corresponding to the 
virial velocity of these massive halos, i.e. about 200 — 300 km s _1 . Since this is typical of the rms velocity of 



6 



a medium sized group, Method 3 turns out to correspond closely with the merging scheme already done in 
NKP94 (although our motivation then was actually to correct for spatial resolution, not overmerging). As 
we will see, the NKP94 results are quite close to the results of our more careful analysis here. We regard 
Methods 2 and 3, and especially Method 1, as less realistic than our adopted procedure. The assumptions 
for each of our breakup methods are summarized in Table 1. 

Table 1. Summary of Overmerger Breakup Methods 



Method 
Preferred 



Method 1 



Method 2 



Method 3 



Assumptions 

M hu = 7 x 10 11 M Q , M inside Sp/p > 1360 
Schechter fragment masses M; M* =< M/L > s ; m L 



CfA> 



*CfA 



-1.26 



v fragment =Gaussian (c,«J m ), <7=rms of galaxies within 1 Mpc or nearest 4 
Fragments positioned by mass; higher mass closer to center 
Brightest fragment mass Mtj = (Mt u M e j j) 1 ^ 2 , at overmerger center, 
and velocity v bf = 0.25 x vj ragment 

M hu = 7 x 10 11 M Q , M inside Sp/p > 1360 

No. fragments from ESD Fig 12a; masses, if ordered, each 90% of previous mass 
v fragment =Gaussian (c,«J m ), <7=rms of 10 nearest galaxies 
Fragment positions not mass-dependent 

M hu = 7 x 10 11 M Q , M inside Sp/p > 1360 

Schechter fragment masses M; M* =< M/L > s ; m L* C j A , a = oiesd = —1-35 
fragment =random Gaussian (<r,wj m ), a = (GMgff/reff) 1 / 2 
Fragment positions not mass-dependent 

Leave simulation overmerges alone 



Merge CfAl grouped galaxies within 



pro] 



235 Kpc , i.e. NKP94 results 



fragment positions enforced 2-cell limit, for all methods 
| v om = velocity of original overmerger 



Figure 2 shows the 1-cell mass distribution of the halos. The no breakup halos follow a power law 
of slope d log N/d log M = —1.33. We stress that the 1-cell and 3-cell masses, being defined by arbitrary 
boundaries, are not particularly physical measures of dark matter halo mass, and in our analysis are used 
only as stepping stones to luminosity. All breakup methods lead to a fairly sharp cutoff in masses at the 
upper end. Note that Method 1 makes an especially noticable pile of fragment masses just below the breakup 
mass limit Mi, u . 

All of these breakup schemes may in fact overestimate the number of fragments. Other purely dissipa- 
tionless simulations (e.g. Carlberg 1994) find that dense DM cores are suprisingly persistent within virialized 
clusters, suggesting that overmerging maybe less significant than generally believed. Also, if M*/10 is equiv- 
alent to a luminosity of L* /10, ESD and Evrard (1994) find too many fragments by an order of magnitude, 
when compared to observations. Our preferred procedure produces even more fragments (a median of 8 
fragments per overmerger for both CDM simulations and 14 for both CHDM simulations, with a maximum 
of ~ 140 for all). Our luminosity assignment method (see §6) prevents overpopulation of visible galaxies, but 
perhaps does not prevent too high a fraction of visible galaxies which are fragments. In any case, as we'll 
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see, our favored statistics turn out to be insensitive to reasonable breakup assumptions, and quite similar to 
our no breakup results. 

4. Properties of Box Groups 

We first describe groups selected using full 3D spatial information from a complete sample of all galaxies 
above a rather low 1-cell mass cutoff. Our procedure for constructing redshift space "sky catalogs" that can 
be compared to CfAl data are described in §6. A group is defined as the galaxies within a bounding surface 
of constant number overdensity. Our grouping scheme is the standard "friends-of-friends" algorithm. The 
critical link distance D is found from the N galaxies in the box above the mass cut by D = Dq(1 3 / N) 1 ! 3 , 
where / is the length of the box and Do is the fraction of the mean interparticle spacing, a dimensionless 
parameter. For each galaxy, all neighbors within a distance D are linked. Each of the newly linked neighbors 
is in turn searched, until no more members are found. We adopted Do = 0.36 (corresponding to a number 
overdensity for selected groups of ~ 20), for our comparisons, yielding a fraction grouped of ~ 70% for CHDM 
and ~ 60% for CDM. This overdensity limit best corresponds with the redshift selected groups described 
here and in earlier work (e.g. Huchra k Geller 1982, Nolthenius k White 1987 (NW), RGH, Nolthenius 1993 
(N93)). In constructing groups near the box boundaries, the periodic images of galaxies were also considered 
during linking. This avoided the problem of clipping groups near the boundaries artificially. 

Table 2 shows properties of groups made from the full box. The CHDM (6p/p) cu t = 80 and CDM 
(6p/p) cu t = 156 boxes will later be assigned luminosities and generate observational sky catalogs, as these 
halo density cuts best match CfAl galaxy densities (see §6). We refer to these below as our fiducial boxes. 
To facilitate comparison between CHDM and CDM at the same halo mass density cut, we also include 
(6p/p) cu t = 80 CDM box results below. Note that these box groups are made from a sample containing a 
high fraction of low mass galaxies and groups, unlike the later groups to be made from magnitude limited 
sky catalogs. 

Table 2. Properties of Box Groups 



simulation 


CHDM 2 


CHDMi 


CDM1.5 


CDM1 


halo (6p/p) cut 


80 


80 


156 


156 


N* 


11305(8585) 


10898(8134) 


14636(8503) 


16050(7247) 


N** 

om 


157 


138 


571 


658 


V c om km s" 1 


309 


316 


315 


322 


V°% gh km s" 1 


283 


319 


266 


291 


N~ grps 


640(620) 


575(575) 


736(680) 


617(514) 


(M v i r /MDM)med 


1.26(1.02) 


1.43(1.18) 


1.27(1.40) 


1.44(1.96) 


/grt 


.71 


.73 


.58 


.64 


DM frac in groups]: 


.15(.23) 


.14(.21) 


.23(.38) 


.27(.50) 



quantities in () are for no-breakup case 

* number of halos in the box 

** number of overmerged halos 

| fraction of halos which are in groups 

I fraction of total DM which is inside the mean harmonic radius of groups 

For our fiducial boxes, the mean spacing between groups is ~ 9 — 12 Mpc for all simulations. However, 
groups are strongly concentrated along filaments surrounded by large voids (see Brodbeck, et al. 1994), so 
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that the distance to the nearest neighboring group is generally smaller; 5 Mpc on average, and 4.7 Mpc 
median. 

Figure 3 shows the average distribution of mass around box groups in all simulations. Figure 3a also 
shows separate densities of the hot and cold fractions for CHDM2. Since virial radii ranged from ~ 0.6 to 
~ 8 Mpc, we minimized smearing in distance by binning the mass about each group as a function of the non- 
dimensional radius r/r^. The result is an unweighted average for all groups. The profile is approximately 
exponential p oc e~ r l rh . Beyond the core, a power law of slope —2 (i.e. an isothermal sphere) fits the 
CHDM curves reasonably well. Simulations of galaxy clusters in a range of cosmologies (Crone, Evrard and 
Richstone 1994) find SI = 1 cosmologies also show an isothermal DM distribution, while £1 ~ 0.2 cosmologies 
do not. CHDM2 densities are ~ 10 — 15% higher than CHDMi. CDM densities, shown in Figure 3b, are also 
exponential. Groups made of halos over (6p/p) cu t = 156 have densities on average ~ 70% higher than those 
made of halos over (6p/p) cu t = 80. The density distribution shows no detectable change of slope beyond the 
radius defined by the halos. Indeed, out at 2r^ the density is still an order of magnitude above the critical 
density. This mass is likely still strongly bound to the group. Figure 4 shows the cumulative mass for the 
fiducial boxes. Note that there is still a significant amount of mass from < r < 2r^ which appears to 
be part of the groups; M(2r h ) is 2.7M(r h ) for CHDM and 2.0M(r h ) for CDM. CDM box groups appear 
to contain a larger fraction of their DM mass at small r. Figure 4b shows the cumulative mass density 
distribution around sky groups. Note that smearing due to grouping in redshift space all but erases the 
differences in the M(r) trends between simulations. Sky groups are on average much brighter and more 
massive than box groups, due to the magnitude limit. These more massive groups have an even broader 
distribution of surrounding dark matter, and M(2rh)/M(rh) is 2.8 — 3.0 for all simulation groups. 

In the Cold + Hot Dark Matter picture, the cold particles fall into gravitational potential wells first, 
later to be followed by the hot particles after they cool. This offset in time continues until the present and 
leads to a higher fraction of cold particles in the cores of the groups, as shown in Figure 5. Interestingly, 
the CHDMi simulation's higher power on large scales seems to lead to a higher contrast between cold and 
hot densities, perhaps through earlier collapse and higher concentration of cold particles today. Data inside 
195 Kpc is unresolved and not plotted. 

We also calculated the virial masses, DM mass content, and velocity bias for these groups. These will 
be discussed in connection with determinations of SI in §8. 

5. Merging the CfAl Catalog 

Several issues need to be considered before we can make realistic comparisons between simulations and 
CfAl data. First is their differing spatial resolution. The galaxy identification scheme requires that a cell 
identified as a galaxy be at a local density maximum. Thus, the nearest possible neighboring halo will be two 
cells away. We require the CfAl data to show the same resolution, on average, before we can reliably identify 
equivalent groups. We therefore attempt to merge CfAl galaxies which are sufficiently close. In NKP94 we 
used a simple merging scheme which assumed an isotropic orientation for the separation vector between all 
galaxies within a group and merged all galaxy pairs separated on the sky by less than 2r c ;/7r, where r c ; is 
the spacing between nearest possible neighbors averaged over the positions of all possible nearest neighbors, 
and 2/tt projects this onto the sky. The isotropic assumption is actually incorrect. Isotropy would only be 
true if the depth of the group were as small as r r j. Here we perform a more careful merging which more 
closely models the resolution of the simulations. As such merging may be important for future comparisons 
between simulations and observations, we describe it in some detail below. 
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Consider a simulation cell tagged to be a galaxy halo. The nearest possible neighboring halo will lie on 
a cube 5 cells on a side ("neighbor cube"), such that inside this cube the 26 cells immediately adjacent to 
the central halo in question cannot be halos. To get an estimate of the mean separation between nearest 
neighboring halos, one could simply average the lengths of the separation vectors between the parent halo 
and all cells on the neighbor cube. To improve this estimate, we first weight each cell on the neighbor cube by 
the frequency that it is actually occupied in the simulations. For example, we find the closest neighbor cells, 
i.e. those offset from the central cell in one coordinate only, are ~ 36% more likely to be occupied than those 
on the edges and corners (average over all simulations). The resulting mean 3D separation between closest 
neighbors is r c ; = 0.509 Mpc . Now consider two concentrations of dark matter within the simulations. If 
the separation r between the centroids of these concentrations is r < r c //2, our halo finder will define these 
as a single, merged halo. If r > r c //2, they will be identified as two halos separated by r r j. Thus, one should 
seek to merge CfAl galaxies closer than r c ;/2 (those between r c ;/2 and r c ; are balanced, on average, by those 
between r c ; and 1.5r c /). 

The last complication is that for real galaxies only sky-projected separations are known. Clearly group 
members in the simulation may be arbitrarily close together on the sky, given the right projection. The 
probability that they are actually close in 3D depends on the depth of the group relative to the pair's 
separation on the sky. Since candidate mergers will have sky separations of order r pro j < 0.25 Mpc and a 
typical group's depth is of order 1.5 — 2 Mpc , the probability of merging is actually low and the isotropic 
estimate of NKP94 results in too many CfAl mergings. To model the probabilities properly, we constructed 
a suite of artificial groups of various sizes (parameterized by the mean pairwise separation of the member 
galaxies r p ), with randomly positioned galaxies. Their density followed roughly the density of halos in 
groups; an r < 0.55 Mpc isodensity core surrounded by a p oc r -2 envelope, truncated so the resulting 
group boundary had an aspect ratio of 1:1:2 with the long axis along the line of sight. This elongation 
corresponds to the observed median aspect ratio of simulation groups in redshift spacce at the fiducial link 
parameters. We used these simulated groups to calculate the probability P me rge(f p , r pro j) that the pair 
should be merged by tallying those pairs close in 3D as well as in sky projection. The resulting probability 
curves vs. r pro j could be well fit with a set of second order polynomials. The output of this exercise was 
then an interpolation table of polynomial coefficients for each of 10 r p 's spanning the range of observed CfAl 
group r p 's. For a median sized group, two galaxies separated on the sky by r < r r j had a probability for 
merging of order 0.1. We then looked at each N93 CfAl group and binary galaxy pair, rolled the dice, and 
either merged or left alone each pair with r pro j < r r j. After several trials, we adopted one which was typical, 
resulting in 39 group members and 5 binaries being merged. Note that this is only 2% of the galaxies. 

6. Constructing Sky Projected Redshift Catalogs 

For each galaxy in the fiducial simulation boxes, radial velocities v r in direction r at distance r were 
calculated from 

v r = (v - v* ) ■ r + H r, (1) 

where the the observer is placed on a "home galaxy" with peculiar velocity vo. For the CHDMi and CHDM2 
simulations, we chose 6 home galaxies satisfying the following conditions: (a) the local galaxy density in 
redshift space (V < 750 km s _1 ) is within a factor of 1.5 of the merged CfAl galaxy density (though still 
usually on the low side), and (b) the closest Virgo-sized cluster is 20 Mpc away in distance. The eight most 
massive clusters ("Virgo" 's) in the CHDMi box were centered about, and within a factor of 2 — 3 of, the 
adopted Virgo mass of M virgo = 2.5 x 10 14 M Q . Corresponding CHDM 2 clusters were less massive. Since 
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the CDM1 and CDM1.5 boxes had the same initial random wave amplitudes as CHDMi, the locations of 
large clusters and filaments were essentially the same as for CHDMi. In order to further remove noise due to 
cosmic variance when comparing the different cosmological models, we chose the 6 viewing locations within 
CDM1 and CDM1.5 to be on the halos nearest to the viewing halo coordinates found for CHDMi. The 
resulting home galaxies had an average peculiar velocity of V pec ~ 800 km s _1 . We retained all data out 
to V = 12000 km s _1 . Since our 100 Mpc box corresponds to 5000 km s _1 for Ho = 50 km s _1 Mpc _1 , 
this required periodically replicating each simulation box on all of its faces. However, at greater distances 
only the brightest galaxies are retained in our magnitude limited catalogs, so the replica boxes are sampled 
sparsely. 

To simulate boundary effects yet still retain a large sample of the sky, we excised a 20° wide "zone of 
avoidance" , leaving a solid angle of A = 10.384 sr. For some comparisons, it is important to match the CfAl 
Survey geometry, despite the much poorer statistics. We therefore generated catalogs which had the local 
"Virgo" rotated to its CfAl right ascension and declination, and the CfAl latitude and declination limits 
imposed, giving a solid angle of A = 2.66 sr. The 10.384 sr catalogs are hereafter referred to as the "sky 
catalogs", the 2.66 sr catalogs are referred to as the "CfA-sky catalogs", and the merged CfAl catalog as 
simply "CfAl". 

An unavoidable problem in comparing all numerical simulations of large scale structure with real ob- 
servations is how to assign luminosities to masses. How nature does this is still understood only in outline, 
and it is not obvious how our 1-cell masses should be assigned optical luminosities. We have therefore made 
the simple but plausible assumption that the baryonic mass in each 195 Kpc 3 cell turns into stars in such 
a way that the resulting luminosity function <f>(L) has a Schechter form similar to that of the merged CfAl 
catalog, and that galaxy luminosity increases monotonically with galaxy mass. 

In solving for the "best" luminosity function of the CfAl catalog, Nolthenius (1993) first assigned 
distances to each galaxy based on the Burstein-Faber (1987) flow model. Others have used a Virgocentric 
infall model. Here, for the purpose of comparisons, it is only important that <f>(L) be defined in an equivalent 
way between the simulations and observations. We therefore assumed simple, unperturbed Hubble flow 
distances D = V/ Ho (except that radial velocities V < 300 km s _1 are set to V = 300 km s _1 ). We next found 
the Schechter parameters a and M* for the merged CfAl catalog using a code based on the inhomogeneity- 
independent method of deLapparent, Geller & Huchra (1989, DGH). <f>* was constrained by requiring that 
the integrated number of galaxies N over the solid angle A of the catalog inside V = 12000 km s _1 matched 
that of the CfAl, 



where Lu m (V) is the luminosity of a galaxy of m = 14.5 at redshift V, and V/Ho is in Mpc. Our <f>* 
constraint differs from the maximum liklihood methods (e.g. Efstathiou, et al. 1990) or that of DGH, yet 
reproduces their unmerged CfAl <^*'s well. After some experimenting, we found that setting m = 14.5 at 
V = 300 km s" 1 for a halo of (6p/p) cut = 80 (CHDM) and 156 (CDM) produced a total density of halos 
~ — 20% above that of the merged CfAl catalog. We chose (bp/ p) cut's of 80 and 156 in order to insure 
that the same small fractions were culled from both the CHDM and CDM catalogs to bring their number 
densities in line with the CfAl catalog. Still, the difference between (6p/p) cu t = 80 and 156 has only a 
minor effect on the sky catalogs since the imposition of an apparent magnitude limit means that only a 
few percent of sky catalog galaxies have 6p/p < 156. (Contrast this to the situation with the box groups, 




(2) 



X(V) = L lim (V)/L* =dex [A(M* - 14.5 + 25 + 5 log(— ))]. 



(3) 



o 
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which are volume limited. The steep mass function then produces a box sample dominated by low mass 
halos and their resultant groups.) Random culling down to CfAl density then yielded our final sky catalogs. 
The 10.384 sr sky catalogs each had 9204 galaxies inside a redshift of 12,000 km s _1 . All 2.66 sr CfA-sky 
catalogs had 2358 galaxies within 12, 000 km s _1 . Because of the periodic boundary conditions, most bright 
galaxies appear more than once in the final sky catalog. On average, each galaxy appears ~ 3 — 4 times in 
the full 10.384 sr sky catalogs, and ~ 1.8 times in the 2.66 sr sky catalogs. This would be a concern for any 
study focusing on large scale structure. However, we focus on median group properties and these should be 
insensitive to this level of replication. 

Our adopted breakup procedure produced a high fraction of fragments in the sky catalogs; 27% (CHDM), 
50% (CDM1.5) and 61% (CDM). Using breakup Method 1, a high percentage of sky catalog galaxies were also 
fragments of overmergers: 27% for CHDM and 70% for CDM. These percentages were only 11% (CHDM), 
37% (CDM1.5) and 46% (CDM) for Method 2, which had many more low mass fragments. 

The "true" mass of a DM halo is not a well defined concept. The DM around a simulation galaxy 
merges smoothly with the DM in a group, which merges smoothly with a general background. We believe 
it's plausible, however, that regardless of how one chooses to define it, the DM mass will have a monotonic 
relation to our 1-cell mass. If we further assume that DM mass and blue luminosity are monotonically 
related, then our luminosity assignments should insensitive to the detailed relations, and be reasonably 
representative. Nevertheless, the relation between DM mass and 1-cell mass is likely to be non-linear, so 
that one would be suprised if a constant M/L resulted from our luminosity prescription, even if in some sense 
M/L is a constant for galaxies. Figure 6a shows M/L for all of our sky catalogs. While it diverges upward at 
the smallest masses (which comprise only a tiny fraction of the sky catalog members), for the great majority 
of galaxies M/L ~ 10 for CHDM and M/L ~ 20 — 25 for CDM. (Other luminosity assignment methods have 
been tried in earlier dissipationless simulations. For example, Gelb (1992) uses the Tully-Fisher and Faber- 
Jackson relations. He too finds a luminosity function which is roughly constant over intermediate masses 
and deviates at the low and high ends.) When the massive halos are broken up, M/L for the fragments ends 
up lower than for the no breakup case. Figures 6(b) and 6(c) show the distribution of 1-cell masses in the 
sky catalogs. Note that breaking up the overmergers changes the mass distribution substantially. A strong 
peak occurs for the fragments, just below the cutoff mass Mj, u = 7 x 10 11 Mq . This is especially true for 
CDM, which has more and higher mass halos than CHDM. 

Table 3 shows the characteristics of the resulting sky catalogs. Note that merging lowers CfAl's lumi- 
nosity density and brightens M* compared to the unmerged CfAl, though the effect is very slight with only 
44 mergers. This is because new mergers which would ordinarily be formed from galaxies with m > 14.5 are 
missing from the CfAl sample. 



Table 3. Sky Catalog Parameters 



Catalog 

CfAl (full) 

CfAl (merged) 

CHDM1 

CHDM2 

CDM1.5 

CDM1 



Y pec 



852±99 
595±247 
700±127 
991±169 



a 
-1.27 
-1.26 
-1.24 
-1.20 
-1.19 
-1.19 



M* 

-21.08 

-21.10 

—21.23 

—21.23 

—21.19 

—21.19 



| peculiar velocity of home galaxy and la range over 6 viewpoints 
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To test the sensitivity of our comparisons to luminosity assignment, we considered two alternate meth- 
ods. Both of these alternates were performed on the no breakup cases only. The first was suggested by J. 
Ostriker (private communication 1993). Cen and Ostriker's (1993; CO) hydrodynamic simulations lead to a 
relation between CDM total density p to t (dark matter + gas/stars) and collapsed baryonic density pbaryon 
(presumably stars). If p = p/p then CO find log(/Sbaryon) = A + B log(/5 to t) + C log 2 (/5 to t)- B and C depend 
mildly on smoothing scale R. In our simulations, we considered p to t the density within a 3-cell cube centered 
on the halo, giving a smoothing scale of R=0.6 Mpc. Extrapolating CO's trends for coefficients B and C 
down to R = 0.6 Mpc then gave B = 2.3, C = —0.20. We assume Mi, aryon / L is a constant. The normal- 
ization A was set to insure that after magnitude limiting (J) p/p) cut = 80 galaxie halos to m = 14.5 for all 
simulations, the number of galaxies retained was similar to the merged CfAl catalog, giving A=6.23 for both 
CHDM models and A=5.8 for both CDM models. This relation, with our distribution of masses, produces 
a Schechter-like luminosity function, but with M* two and a half magnitudes too bright (~ —23.5), and a 
too steep (~ —1.7). The resulting halo M/L (see Figure 6a), attaches much higher L's to massive galaxies, 
relative to the Schechter prescription. This leads to a disproportionate number of distant, luminous galaxies 
at the expense of those closer. 

The second method uses the blue Tully-Fisher relation from Fouque, et al. (1990). We assume the cell 
mass M is gravitationally bound so that a circular velocity may be defined 



GM 

Vcirc = ] J — . (4) 

The Fouque, et al. Tully-Fisher relation then defines the absolute B magnitude Mb of the galaxy as 

M B = -5.5log(y eire ) - P (5) 

where (3 is a calibration. This is equivalent to M 11 / L = const. Using the Fouque, et al. value of (3 = 8.0 
corrected by —0.37 for average internal extinction in spirals and by —0.29 to put on the Zwicky system results 
in over an order of magnitude too few galaxies surviving the m = 14.5 limit! This isn't too suprising. Our 
halos are 200 Kpc across, ~ 10 times bigger than the visible size of a typical spiral, and the relation between 
our 1-cell mass M and the observational M(r ~ 10 Kpc ) is not known. We therefore adopt a [3 which best 
matches the merged CfAl galaxy density of 886 sr _1 . The resulting halo M/L attaches much higher L's 
to low mass galaxies, relative to the Schechter prescription. This leads to catalogs with a disproportionate 
number of nearby galaxies and strongly favors fainter, lower v gT groups. 

7. Tuning the Grouping Algorithm Link Parameters 



Conventionally, a group is thought of as a gravitationally bound collection of perhaps 3 to 50 galaxies 
which has collapsed. For our purposes, we consider a group more generally as any set of galaxies which 
satisfy a linking criteria on the sky and in redshift. Any valid cosmological model must be able to produce 
a set of such groups whose typical properties are in good agreement with identically selected real groups for 
any link lengths. It is not necessary that the groups be bound, let alone collapsed or virialized. Nevertheless, 
it is also of interest to identify that set of groups which is closest to satisfying the conventional definition 
with the least contamination by interlopers, e.g. for estimating £1 from the M/L method (see §8). The 
corresponding link criteria which make these most realistic groups are referred to below as the "fiducial" 
links. 

We use the adaptive grouping algorithm of Nolthenius (1993), which is a modification of the original 
friends-of-friends algorithm in redshift space given by Huchra and Geller (1982). Galaxies at redshift V 
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km s 1 are linked if the separation on the sky is less than Dl and their difference in redshift is less than Vl, 
where 

< f uW = Vr(i + «,W), (6) 

7T 



D L (V) = D n R-y s <V )^) , (7) 

A(V) is given by equation (3). Vo is an arbitrary redshift distance used for scaling, here taken to be 
1000 km s _1 . D n , an input parameter, is the sky link expressed as a fraction of the sky-projected mean in- 
terparticle spacing Ft^i s ■ ^(^0 ^ s ^ ne integrated galaxy luminosity function above the apparent magnitude 
limit visible at redshift V. Dl is designed to insure that a group at the minimum number density contrast 
just meets the selection criteria at all redshifts (see NW). We scale Vl with the mean interparticle spacing. 
We chose this scaling, rather than the linear NW scaling (whose resulting groups j^(V) trend turned out 
to match slightly better to groups selected using full 3D information at the same overdensity limit), simply 
because it produced a flatter distribution of group M/L with distance. We have also run cases with the NW 
version of the algorithm and find negligible differences in group properties and no change in our conclusions. 
For the purposes of comparing models, either V(L) can be used. We parameterize the size of the link in 
redshift with V 5 = Ft (5000 km s" 1 ). 

Consider a collection of galaxies which are close together in redshift space. Their rms velocities about 
their mean peculiar velocity will govern how elongated along the line of sight this group appears, providing 
a diagnostic for determining the optimal redshift link V5. Define the aspect ratio A z of this "finger of God" 
seen in redshift space as A z = V(ARA + ADec)/2AV , where A refers to the maximum extent of the group 
in right ascension, declination, and line of sight velocity. Define a group's aspect ratio A r in real space as 
A r = V(ARA + ADec)/2A(V — V pec ). Note that A r can be defined only for simulation groups, since the 
galaxies' true distances are then known. Because groups are often poorly isolated from their neighbors and 
because of random peculiar velocities, the true depth of groups picked out in redshift space will inevitably 
be larger than the sky-projected group size; groups will be elongated along the line of sight. This is not the 
familiar "finger of God" one sees in redshift maps, it is its counterpart seen in real space. This elongation 
would, on average, disappear if one could assign memberships with perfect knowledge of galaxies' true 
distances. Figure 7 shows the trend of A r vs. V5 for our simulations, at the fiducial D n = 0.36. Consider 
a collection of galaxies all close together in redshift space and with a small rms velocity v gT about their 
mean peculiar velocity. At V5 below some threshold, one picks out only groups whose members are within 
~ Dl of each other in real space, and A r is constant and near a minimum (in fact, since there are two sky 
dimensions and only one depth dimension, there is a tendency for A r to rise slightly at very low V5, when 
only a fraction of the valid members are being selected). At V5 below the threshold too many valid group 
members are excluded, and we refer to this as the "clipped regime". Above the threshold V5, A r begins 
to rise as outliers begin to significantly contaminate the group. High V5 is then said to produce groups in 
the "interloper regime". For the CHDM groups, this occurs at V5 ~ 300 km s _1 . The CDM curves never 
flatten, showing that groups have significant contamination even at the lowest V5. The CDM1.5 curve slope 
steepens slightly at V5 ~ 500 km s _1 when the rate of contamination rises. The CDM1 curve shows only a 
very gradual curvature and a slight discontinuity in slope at V5 ~ 800 km s _1 . These curves suggest fiducial 
V5 ~ 300 km s _1 for CHDM and ~ 600 — 700 km s _1 for CDM. Note that the minimum A r does not approach 
1, due to the inherent loss of depth information in redshift space. It is possible to force rounder groups and 
thus smaller j4 r 's by raising D n and restricting V5, but at the cost of rejecting an unacceptably high number 
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of valid members. For comparison, Figure 7 also shows the curves for groups from a sky catalog made from 
a 100 Mpc box of Poisson distributed particles with Gaussian random velocities of a = 200 km s _1 and 
a = 350 km s _1 . In these catalogs, there is no coherent clustering to confine the depth dimension, which 
therefore rises steeply with V5 even at small scales. Note by comparison that even the CDM catalogs show 
significant coherent motion by their more slowly rising curves. 

Perhaps the most relevant measure is to ask what V5 is required in order to produce the same median 
i> gr as is seen in groups selected using full 3D information on our fully corrected simulations at the same 
fiducial D n = 0.36. This is shown in Table 4. This is V 5 ~ 400 km s" 1 for CHDM, V 5 ~ 550 km s" 1 
for CDM1.5, and 700 km s _1 for CDM1. If no breakup is performed, i> gr 's are very similar, but require 
somewhat lower Vs's to match up with 3D groups. 

Table 4. V 5 Link Best Matching 3D Groups' v gT V 

simulation CHDM 2 CHDMi CDM1.5 CDM1 

V 5 398(371) 413(330) 551(456) 702(575) 

(vgr ) med 126(151) 130(148) 187(188) 241(237) 

*Vs link giving redshift space groups with the same (i> gr ) me( j as 3D groups, 
at fiducial D n = 0.36. All in km s _1 , no breakup case is in parentheses 



Finally, we note that the fraction of galaxies grouped / gr will at first rise steeply with V5 as new 
members are rapidly added from the dense region containing the group. The / gr (V5) curve will then show 
a distinct drop in slope as it enters the interloper regime, when essentially all valid members have been 
added and primarily outliers and interlopers from the low density surroundings are incorporated at higher 
V5. This transition occurs near V5 ~ 350 km s _1 for both CHDM simulations, and near V5 ~ 600 km s _1 for 
both CDM simulations (see §9). As these values are representative of the other measures above, we adopt 
V 5 = 350 km s" 1 as the fiducial redshift link for CHDM, and V 5 = 600 km s" 1 as that for CDM. 



8. The Dynamical State of the Groups and Implied £1 



To assess the dynamical state of groups, we calculated the number- weighted virial mass estimator for 
3D-selected box groups. This is the M„; r appropriate if the mass is dominated by a dark matter background 
distributed like the galaxies, and is given by 



M vlr = H^L, (9) 

where, for a group of n galaxies, a = (n — 1) _1 — * s ^ ne 3D rms velocity about the mean velocity 

v of the group. We then compared this to the true mass Mom in dark particles within the mean harmonic 
radius r^, where 



n(n — 1) . 
r h = I j_ , (10) 

and rij is the separation between galaxies i and j. 

Figures 8a and 8b show scatter plots of group log M v i r vs. logMpM for the no breakup and breakup 
CDM boxes, respectively. Many of the no breakup CDM groups appear to be unbound, with the median 
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M v i r / Mum as high as 2.0 for CDM1, while the breakup box appears only slightly biased to high M v i r . This 
effect is not due to the higher minimum halo overdensity ((6p/p) C ut = 156 vs. (6p/p) cu t = 80 for CHDM), 
since it is also clearly seen in the (6p/p) cu t = 80 boxes. It also appears equally strong when only better 
sampled groups (at least 10 members) are considered. The explanation is that many small groups do not 
appear in the no breakup cases, because of overmerging. Recall that a halo is defined as a local density 
maximum, requiring neighboring cells to be relatively underdense. Also, a group must contain at least 3 
members to be identified. The overmerged CDM halos often dominate their surroundings to such an extent 
that too few other local density maxima can be defined, preventing a group from being identified. When 
these are broken up, new groups result. These groups turn out to have M„; r close to Mom, on average, so 
that their inclusion lowers the median M„; r / 'Mom nearer to unity. This issue affects CDM much more than 
CHDM, since CDM has ~ 4 times as many overmergers (see Table 2). 

Note that the breakup case in Figure 8b still has these high M„; r / 'Mom groups. Correlating high 
M v i r I 'Mom in the no breakup CDM box with other group properties, we found that these groups tend to be 
poorly sampled (all groups with M„; r / 'Mom > 30 have 5 or fewer members), have higher than average size, 
low density contrasts, high peculiar velocities, and tend to lie in denser regions. Using computer visualization 
and color-coding groups by their M„; r / 'Mom shows these groups tend to be small and lie near the outer 
parts of dense clumps and filaments. Putting this all together suggests these groups may have experienced 
tidal shearing sufficient to unbind them. Other recent work supports this idea (e.g. Mamon 1994). Some 
high Myir/MoM^ are no doubt simply due to the inherent high noise in M v i r . Figure 9c shows M„; r vs. 
Mom for groups from all 6 fiducial linked sky catalogs. Note that the high V5 required to match median 
3D i> gr 's nevertheless produces a tail of spurious groups with high v gT , little mass within r^, and hence 
high M v i r I 'Mom ■ The high V5 link essentially produces a floor for M v i r , which is unreasonably high for the 
smallest groups. This is a general feature of all group catalogs linked in redshift space (see Geller & Huchra 
1983; HG, RGH). Banding is due to using only a 10% sample of the dark particles; the smallest (spurious) 
groups contain only a few dark particles. 

The corresponding plots for the CHDM boxes are in Figure 9. CHDM box groups appear little changed 
by breakup, partly because they contain only a quarter as many overmergers, and partly because they show 
lower random velocities on small scales. Note in Figure 9c that the lower V5 level appropriate for CHDM 
significantly reduces the problem of spurious low Mom but high M„; r / 'Mom groups. Since a low V5 also 
appears appropriate to the real universe, the problem of contamination and spurious groups may not be as 
severe as seen e.g. in Geller and Huchra (1983) 

Note that massive groups (M > 10 14 Mq ) in all boxes appear to define a narrow M v i r / Mom which 
appears clearly virialized, yet which lies systematically higher than the virialization line. There are two 
biases towards high M„; r / 'Mom which likely account for this, and which are present for all groups. First, 
the r/j sphere within which particles are counted is centered on the mean position of the halos. Centering 
instead on the dark particle concentrations lowers M„; r / 'Mom ■ We confirmed this by centering on the mass- 
weighted mean position of the particle condensations within 1.5r/j of the groups and found it raises Mom 
such that the median M„; r / 'Mom is lowered by 7%. Second, many groups, especially the larger ones, are 
elongated density enhancements along filaments. Counting within a sphere will include spurious low density 
areas along the "equator" while missing the high density regions along the "poles" . If this bias is ~ 30% 
in Mom, which seems plausible (see group geometries in the video sequence of Brodbeck, et al. 1994), then 
M > 10 14 Mq groups in the breakup boxes appear to obey the virial theorem, on average, in all simulations. 

We also determined several measures of the "velocity bias" b v (Carlberg & Couchman 1989): the velocity 

of the galaxies compared to that of the underlying cold dark matter particles. The global velocity bias b 9 v lohal 

1/2 

is defined by a sum over the i galaxies and j DM particles within the box b 9 v lobal = (J2i v f ga \l v ] dm) 
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The results are summarized in Table 5 below and shown in Figure 10. These change significantly when one 
looks only within the virial radii of the individual groups. The velocity bias within a group's is defined 

1/2 

over the i galaxies and j DM particles as h a J v = (J2i( v i ~ Vgai) 7 ' /J2j( v j ~ vdm) 2 ) , where v ga \ is the mean 
unweighted galaxy velocity and vdm is the mean DM particle velocity within the group. CHDM galaxies 
show a stronger bias within groups. CDM1 galaxies, however, show little or no bias, in or out of groups, 
in agreement with b 9 v lohal = 0.94 from Katz, Hernquist, & Weinberg's (1992) hydrodynamic simulation. 
CDM1.5 shows a moderate bias both globally and inside groups. 

Table 5. Velocity Bias 



simulation ¥ v lohal < WJP > avg < WJP > med 
CHDMi 0.93(0.90) 0.76(0.78) 0.78(0.76) 
CHDM 2 0.92(0.89) 0.72(0.69) 0.72(0.71) 
CDM1.5 0.83(0.79) 0.76(0.74) 0.75(0.79) 
CDM1 0.95(0.90) 0.97(0.92) 0.98(1.06) 
No breakup case is in parentheses 



The standard method for estimating £1 from bound virialized groups is to assume all galaxies have the 
same mass-to-light ratio M/L, given by the median M/L for groups, then integrating over the luminosity 
function to get the mass density (Kirschner, Oemler & Schechter 1979). For a Schechter luminosity function 
this means 

" = ||<f )4>*L*T( a + 2). (11) 

Table 6 shows the resulting inferred £1 using the fiducial grouping algorithm parameters and Schechter 
L(M) for the sky catalogs. Means and standard deviations are over the six viewpoints. Using the Cen- 
Ostriker L(M) gives very similar Si's. No N(z) corrections (see §9) were applied to get < M/L > me d since 
group M/L is virtually independent of redshift. 

Table 6. SI from Group M/L Method 



catalog V5 

CfAl (full)* 350 

CfAl (merged) 350 

CHDMi 350 

CHDM 2 350 

CDM1.5 350 

CDM1 350 



< M/L > med 
62 
79 

118 ± 19(127 ± 9) 
139 ± 12(120 ± 7) 
159 ±9(164 ±7) 
177 ± 9(186 ± 20) 



SI 
0.06 
0.08 

0.12 ± .014(0.12 ± .003) 
0.12 ± .010(0.11 ± .003) 
0.15 ± .007(0.16 ± .009) 
0.18 ± .007(0.19 ± .012) 



CfAl (full)* 600 119 0.12 

CfAl (merged) 600 152 0.15 

CHDMi 600 198 ± 14(193 ± 10) 0.21 ± .008(0.18 ± .008) 

CHDM 2 600 200 ± 14(178 ± 10) 0.18 ± .010(0.16 ± .007) 

CDM1.5 600 350 ± 23(308 ± 16) 0.34 ± .030(0.31 ± .016) 

CDM1 600 365 ± 24(372 ± 26) 0.37 ±.012(0. 38 ±.022) 

* full, unmerged CfAl catalog. 

* No-breakup results in parentheses 
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Note that all of our SI = 1 simulations yield low observed Si's. Thus, these simulations show a feature 
which has been seen in real data for over a decade and which has been used by some to argue for a low 
SI universe: How can SI be 1 if bound groups and clusters consistently account for only ~ 10 — 20% of the 
mass? Three factors can explain this discrepancy. (1) Only the mass within the mean harmonic radius 
of the groups is measured by the virial estimator. We saw in Figure 3 that density falls approximately like 
an isothermal sphere outside group cores, and that this continues to at least 2r^ (about midway to the next 
nearest group, on average), where the cumulative mass M(2r^) = 3M(r/,) (see Figure 4b). This additional 
factor of 3 in mass is almost certainly bound to the groups, since p(2rh) — 4/9 c (CHDM) to 10/9 C (CDM), and 
infall is still occurring. The fraction of the total mass in the box which lies within the sphere of fiducial 
linked groups is only ~ 15 — 25%, as shown in Table 7. Most mass lies outside of groups. 

Table 7. Mass, Volume Fraction Within Box Groups 



simulation 


N* 

grps 


J mass 


f f 

J v olume 


CHDMi 


575 


0.14 


0.0086 


CHDM 2 


640 


0.15 


0.0072 


CDM1.5 


736 


0.23 


0.0038 


CDM1 


617 


0.27 


0.0026 



* Number of groups in box 

** fraction of box mass which is within of groups 
| fraction of box volume which is within of groups 

(2) The relevant virial velocities appropriate for measuring the local mass are actually those of the individual 
particles, not the galaxies. As already seen in Table 5 and Figure 10, the median velocity bias within CHDM 
groups is ~ 0.75. Since M true /M v i r oc l/b^, this contributes another factor of 1.8 to the SI estimate. (3) 
We've also assumed spherical symmetry in Figures 3 and 4. The true shape of groups is more elongated, as 
structure is still quite stringy at this point in the evolution. The inappropriately counted matter along the 
"equator" of the groups is of lower density and fails to compensate for the uncounted high density regions 
missing above the "poles" by perhaps a factor of 1.3. The product of these factors (~ 7) approaches an 
order of magnitude. While factor (1) above is closely related to the fact that the galaxies are a biased 
(b=1.5) tracer of the DM, it seems likely that the less clustered hot DM, especially that filling the voids (see 
Brodbeck et al. 1994) will additionally bias M/L to the low side. Together, these factors show that, within 
the assumptions of this method, an observed SI ~ 0.1 can indeed be consistent with a true SI = 1 universe. 
In fact, virial estimates appear to be underestimates of the true mass in real clusters as well. At present, 
gravitational lensing appears to be the most reliable method of measuring total masses, and several studies 
show that typically the resulting total masses of clusters out to just beyond the optical radii is roughly a 
factor of 2 or 3 higher than virial estimates (e.g. Kaiser, et al. 1994). There are other, more subtle problems 
with the M/L method; for example, the stellar populations of groups and especially clusters is well known 
to be older and have lower M/L than is typical for the field galaxy populations. 

Our SI at V5 = 350 for the CHDM simulations is 50% higher than that for the merged CfAl at the same 
V5: Qchdm = 0.12 vs. Q,cjai = 0.08. The difference is due almost entirely to the difference in the median 
of the mean harmonic radii of groups, as discussed in §10. Our SI for the full, unmerged CfAl is smaller 
than that of HG's SI ~ 0.1 for the nearby, all sky m=13.2 CfA sample, and smaller than RGH's SI ~ 0.13. 
As we've argued elsewhere (NW, N93), we prefer a significantly smaller redshift link than these authors, 
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which directly lowers the virial masses. The RGH work is for the CfA Slices, whose group's show a suprising 
60% higher v gT than for their equivalently selected CfAl groups, due to differing sample depth and, to some 
extent, cosmic variance (see RGH for discussion). Also, their link parameters differ from ours. They used a 
smaller sky link (giving smaller r^) but larger redshift link, leading to higher v gT 's and net higher M/ Us. 

When using full 3D information to make sky catalog groups (D n =0.36 as before), our median Si's are 
higher: ft =0.17(both CHDM), 0.36(CDM1.5), and 0.48(CDM1). Note that the V^s necessary to match 
3D v gT 's are, except for CDM1.5, higher than our fiducial V5, as shown in Table 4. For example, raising 
CHDM's V5 from the fiducial 350 to 398 km s _1 raises < v gT > me d by ~ 15% and hence < M/L > me d and 
ft by ~ 32%, lessening the difference between the observationally inferred and true ft. The same holds for 
CDM1. Our CDM1.5 estimate of ft ~ 0.33 agrees well with that of Katz, Hernquist & Weinberg's (1992) 
TREESPH hydrodynamic simulation result of ft ~ 0.31 on a much smaller sample. 

9. Sky Catalog Results and Comparisons with CfAl 

When comparing our sky catalog groups with those of CfAl, there is one more calibration to consider. 
To maximize the amount of precious observational data used, we retain the full magnitude limited CfAl 
catalog rather than a volume limited subset. In a magnitude-limited catalog, both the sky and redshift 
"friends-of-friends" links scale up with distance, and group properties will change significantly with distance. 
In particular, more distant groups will be larger, brighter, and have higher v gT 's. Median values will therefore 
be biased if groups are distributed differently with distance between the two datasets. Chance differences in 
large scale structure will mean, in general, that groups are in fact distributed differently with redshift. Our 
home galaxy selection criteria do not fully insure that simulation groups are distributed in redshift like CfAl 
groups. Equally important is our small box size. Periodic boundary conditions give repeating structures 
every 100 Mpc = 5000 km s _1 , whereas the CfAl data are actually rather sparse beyond the Virgo Cluster. 
Figure 11 shows the galaxy density vs. redshift for the merged CfAl and CHDM2 sky catalogs. Relative to 
the CfAl dataset, CHDM2 remains underpopulated out to ~ 3000 km s _1 and overpopulated beyond. The 
other simulations follow this same pattern. Not correcting for this "N(z) bias" will lead to overestimating 
the average or median sizes and v gT 's of simulation groups, as well as affect other properties. We've handled 
this by averaging results for four random subsamples of the simulation groups such that, when their redshifts 
are binned to 1000 km s _1 bins, the number of groups vs. redshift N(z) matches that of the CfAl groups 
selected at the same links. The median group properties we present below which are medians from these 
subsamples are labelled "N(z) corrected" for clarity. 

In Figure 12 we show v gT vs. V5 at our fiducial D n = 0.36. It is analagous to Figure 2 in NKP94, 
which was done at D n = 0.47 to show good agreement in CDM1.5 velocities between our work and that 
of Moore, Frenk and White (1993; MFW). They used a quite different particle-particle-particle-mesh code 
with better spatial resolution but poorer mass resolution than ours here. This figure emphasizes again one 
of our main conclusions; CDM groups have much higher internal rms velocities than observed, while CHDM 
is in good agreement. It may be suprising to see that the curve for CDM1.5 differs so little from that of 
CDM1, since numerous studies show biased CDM has lower velocities. Our global rms peculiar velocities (i.e. 
for all galaxies within the box) are indeed much higher for CDM1 than those for CDM1.5; 944 km s _1 vs. 
650 km s _1 , or 45% higher. Within groups, however, two different effects combine to reduce this difference. 
First, on small scales (~ 1 Mpc), CDM1 rms peculiar velocities are only 30% higher than for CDM1.5. This 
is true both for the full set of box groups, and for the 3D selected groups from the sky catalogs. Second, 
in redshift space, the size of the V5 link strongly affects v gT . Using the same V5 = 600 km s _1 on both 
simulations will then further reduce their differences to only 7% (i> gr = 216 km s _1 vs. 202 km s _1 ) This 
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same effect is seen in MFW, and our result agrees well with the 8% difference seen between their b=1.6 and 
b=2.0 CDM median v gT 's. 

It is important to note that the error bars on these and later curves are la deviations from the 6 different 
sky catalogs. While our intent was to measure cosmic variance, in fact these are a significant underestimate. 
Because of the small size of our box, many of the same groups are seen in most or all viewpoints, (albeit 
sampled differently due to the magnitude limit). The difference between the CHDMi and CHDM2 curves 
is a fairer estimate of cosmic variance. Our estimated probability of a 100 Mpc 3 box having the power 
spectrum of CHDMi of ~ 10% (NKP) corresponds to about 1.7 a. If so, then the difference between the 
CHDMi and CHDM2 curves may be a rough estimate of 1.7 a error bars due to cosmic variance. This is the 
method assumed in quantifying comparisons with CfAl below. 

Another important point is that since CHDMi, CDM1, and CDM1.5 all have higher power in the 
longest waves within the box, their "Virgo" 's are larger and richer, and the fraction of galaxies in groups is 
artificially higher than would be typical. If one believed the CfAl data were a "fair sample" (but see below), 
then comparing to CHDM2 would be more appropriate. To estimate this effect on the CDM curves one can 
do a rough calibration by looking at the CHDM2 curves shown here and below and shift the CHDMi curve 
towards it, carrying along rigidly the two CDM curves in the same direction. In fact, however, the luminosity 
density of the CfAl sample appears to be ~ 25% higher than for the much larger (fair sample?) APM data 
(Tully, private communication). If the presense of rich structures like Virgo, Coma, and the Great Wall in 
this sample similarly indicates slightly unusual higher power on larger scales, then the appropriate curve to 
compare to CfAl may be intermediate between CHDMi and CHDM2. 

Figure 13 shows f gT , the fraction of galaxies grouped, vs. V5. Fraction grouped is a powerful statistic, 
since by separating our grouping algorithm into velocity and sky projected components, f gT is sensitive 
not only to small scale pairwise velocity differences, but also the degree of spatial concentration of galaxies. 
CDM differs from CHDM significantly on both measures. As remarked earlier, the improved merging scheme 
applied to the CfAl data here had the effect of performing fewer mergers. This left more CfAl members in 
dense regions and thus raised the fraction of galaxies grouped; e.g. from 42% to 49% at the fiducial links. 
It had very little effect on the v gT . Another change from NKP94 was to make CDM sky catalogs using 
only halos above (6p/p) cu t = 156 rather than 80. This had the effect of slightly raising CDM's v gT 's by 
~ 3%. The net result is closer agreement between CHDM2 and CfAl f gT curves. CHDMi 's higher power on 
large scales leads to a significantly higher correlation and a too high / gr . CDM's high small scale pairwise 
velocities and puffier, less concentrated galaxy filaments inhibit grouping, especially near the fiducial links. 

By combining Figures 12 and 13 to eliminate the grouping link V5, we can both enhance the differences 
between CDM and CHDM and produce a statistic which is very robust. This, our favored statistic, v gT vs. 
/ gr , is shown in Figure 14. CHDM2 is in close agreement with the CfAl data. CHDMi groups too high a 
fraction of galaxies, while CDM both groups too few galaxies and produces v gT 's too high, at the ~ 4<r level 
(due to cosmic variance) for CDM1 and ~ 5<r level for CDM1.5. The suprisingly high discrimination shown 
by this statistic shows the power of grouping in redshift space. The "pressure" provided by higher small 
scale pairwise velocities will tend to expand spatial structures on ~Mpc (i.e. galaxy group) length scales. 
This not only lowers the fraction of galaxies grouped, but also raises their rms velocities, and does so at all 
link criteria. f gT is sensitive to the presense of large clusters in relatively small samples like the CfAl data. 
It is therefore important to see how v gT vs. f gT behaves on simulation data with the same sky coverage 
as CfAl, shown in Figure 14(b). The dominance of the local "Virgo" 's tends to raise f gT for the CHDMi, 
CDM1, and CDM1.5 simulations (which, recall, have not only higher power, but identical locations for large 
scale structure). CHDM2, which has smaller "Virgo" 's, actually shows a small decrease in f gT . Median 
v gT , as expected, shows little change. The net conclusion remains the same, albeit noisier and with higher 
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cosmic variance. CHDM2 again shows excellent agreement with CfAl throughout the V5 range. Comparing 
data point by point, we see that CHDMi has v gT 's slightly too high, while CHDMVs agree closely with 
observations. CHDMi now seems clearly to group too many galaxies. For a curve intermediate between 
CHDMi and CHDM2 (see §10) Figure 14b can then be interpreted as indicating that group analysis favors 
a slightly lower £1„ than the 0.30 used here. 

The / gr vs. v gT statistic is also quite robust. Figure 15 shows v gT vs. / gr for the no breakup simulation 
sky catalogs, without N(z) correction, for our three L(M) assignment methods. Using the (b) Cen-Ostriker 
L(M) prescription leaves the curves almost unchanged, even though the luminosity function's M* is now 
~ 2.5 magnitudes too bright and the faint end slope is a ~ —1.8. The (c) Tully-Fisher (TF) prescription's 
L(M) differs even more drastically from observations, and is quite un-Schechter-like. Still, the resulting v gT 
vs. / gr curves are qualitatively similar to the Schechter and Cen-Ostriker results, and again CDM curves 
are too high while CHDM's curves are now too low. Note that all simulation TF v gT 's are lower. This 
is because the TF L(M) strongly picks out nearby, low-mass galaxies and groups and therefore gives lower 
median v gT 's. Comparing Figures 14a and 15a shows that breaking up the halos and correcting for differing 
redshift distribution lowers v gT , but still leaves CDM much too high and CHDM2 in good agreement with 
observations. Notice also that the corrected CDM1.5's curve actually lies above the CDM1 curve. The 
reason is that CDM1.5 groups a significantly lower fraction than CDM1. Thus it is more accurate to say 
that the CDM1.5 curve is left of the CDM1 curve. 

Figure 16 compares the breakup methods on the v gT vs. f gT plane. Relatively little change is seen 
for any method when applied to CHDM, as shown in Figure 16a for CHDM2. For CDM1, shown in Figure 
16b, all methods are close except for breakup Method 1. By forcing essentially equal, maximum possible 
masses for all DM halo fragments, Method 1 raises f gT dramatically while lowering v gT . Figures 16a and 
16b Method 1 curves lie virtually on top of one another. In fact, breakup Method 1 has the unfortunate 
property that all simulation curves are almost degenerate on this plane, and significantly below the CfAl 
curve. While there is no reason to think that this breakup scheme is reasonable, it does show that it is 
possible to construct sky catalogs whose v gT vs. f gT properties are not discriminated. 

Figures 17-19 show how v gT vs. f gT change when V5 is held fixed and D n is varied instead. Here, groups 
grow primarily on the sky and only secondarily in redshift depth. In NKP we described Figure 17, which does 
not have the breakup or N(z) corrections included. Our conclusion there was that the dense cores picked out 
at small D n were significantly "cooler" in CfAl than for any simulation. Figure 18 shows that this is only 
true when overmergers are not broken up and no N(z) correction is made for differing redshift distributions. 
When corrected in this way, CHDM actually reproduces the CfAl results very well, while CDM v gT 's still 
remain too high. The largest effect is the N(z) correction. As described earlier, the simulation groups tend 
to lie at higher distance, where the magnitude limit then identifies larger, higher v gT groups. Taking the 
median of random sub-samples whose distribution in distance matches that of the CfAl lowers median v gT 's. 
Comparing Figures 17 and 18 shows that f gT in fact changed very little at these fiducial Vs's, while most of 
the change is in v gT . Figure 19 is similar to Figure 18, but for the CfA-sky catalogs. Again, at these low 
and moderate Vs's, f gT changes very little between these samples. 

Figures 20 and 21 shows the number of groups per steradian vs. redshift link, and vs. sky link, and 
demonstrate the percolation properties of the catalogs. As V5 is raised, new groups will be identified, and 
some existing groups which are close in redshift will be merged. Beyond V 5 = 350 km s _1 , the CHDMi curve 
actually shows a significant decline as merging outpaces the production of new groups. An intermediate curve 
between the two CHDM curves would appear to agree much better with CfAl groups, but appears likely to 
peak at both D n and V5 too low. Comparing such a CHDM curve to the CDM curves, which roughly peak 
together with CfAl, indicates that a lower £1„ would improve CHDM's agreement. Pure CDM, however, 
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appears to significantly underproduce groups at small V5. Thus, on this measure, CHDM appears to agree 
well if CfAl has a moderately high amount of large scale power, and if we lower £1„ (but not to zero). 



10. Possible Problems 



One quantity which shows significant disagreement with observations for all simulations is group size. 
Figure 22 shows the median vs. D n . All simulations are 50% — 80% higher than CfAl, at least near 
the fiducial link of D n = 0.36. CDM actually fits slightly better. This is because CDM groups galaxies 
less efficiently and groups tend to be poorer and smaller than CfAl groups. We've been careful to filter the 
observations to the same spatial resolution as the simulations, so the difference may be real. Or, it may 
instead reflect fundamental limitations in our simulations. One possibility is that the overdense DM cells 
we call galaxies still retain some of the distribution properties of their parent DM, and are insufficiently 
"galaxy- like". Serna, et al. (1994) show that cluster galaxies residing in a DM background will show a 
stronger concentration than the DM, by a factor of about rh(DM)/rh(gal) ~ 2, higher for older clusters. 
If CfAl galaxies are like Serna's galaxies, while our DM cells behave somewhere between galaxies and DM 
particles, one might expect rh(gal)/rh(DM) values similar to what we see. 

Interestingly, our (disfavored) breakup Method 1 gave median values which were closer to that for 
CfAl. This is because such a large number of fragments met the magnitude limit, since we artificially 
imposed nearly equal masses on all fragments. This produces many more very small groups. Even so, the 
trend of median with D n was still steeply negative, similar to Figure 22. Recall that with Method 1 all 
simulations produced virtually the same i> gr vs. / gr curves and in poor agreement with CfAl. 

11. Discussion and Conclusions 



We've shown that COBE-normalized CHDM at £1„ = 0.3 produces group properties similar to those of 
the CfAl, while CDM at b=1.0 and b=1.5 groups too few galaxies and gives i> gr too high, at the several 
a level. We now attempt to refine our estimate for an £1„ which is in optimal agreement with our group 
analysis. 

The CHDMi power spectrum in real space P r (k) is a factor of 2 higher than that for CHDM2 on scales 
comparable to the size of the box (KNP). The redshift space power P z (k) is amplified by a factor 



2 Q° 6 1 Q° 6 

P z{ k) = P r (k)[l + -— + -(— f] (12) 

(Kaiser 1987). For our SI = 1 b=1.5 CHDM simulations then, P z (k) = l.bP r (k), so that in redshift space 
CHDMi has P z (k) 3 times higher than that of CHDM2. The CfAl P z (k) is approximately a factor of 1.6 
higher than that of the CfA2 at these scales (Vogeley, et al. 1992), while it may be more comparable to CfA2 
on smaller scales. The CfA2 P z (k) is in turn a factor of 2 higher than that of the much larger APM Survey 
on similar scales (Baugh & Efstathiou 1993). While these scales are larger than typical group/clustering 
lengths, our CHDMi vs. CHDM2 simulations show that higher power on these scales indeed "crosstalks" into 
higher fractions of galaxies grouped (see Figure 13), perhaps aided by percolation along extended filaments. 
Such coupling of large to small scales has already been noted for pairwise velocities (Gelb, et al. 1993). Also, 
on galaxy scales the CfAl galaxy luminosity density is a factor of 1.25 higher than that of the APM. If the 
APM can be taken as approaching a fair sample of the universe, the CfAl then appears to have a redshift 
space power spectrum at large scales which is at least twice as high as is typical. This then suggests that 
the proper curve to compare to the CfAl data is intermediate between those for CHDMi and CHDM2, but 
closer to CHDMi. Figures 13 and 14 then indicate an optimum £1„ somewhat lower than 0.3. While the 
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detailed relation between £1„ and v gT vs. / gr is as yet unexplored, if it is approximately linear then Figures 
13 and 14 suggest an optimum £1„ ~ 0.2. With the intermediate CHDM curve described above, Figures 20 
and 21 would be in excellent agreement with the CfAl data for any £1„ significantly less than 0.3. Figures 18 
and 19 favor £1„ ~ 0.3, but are, within the errors, compatible with a slightly lower value. The total mass m„ 
of all massive neutrino species is related to £1„ by m v = 23.51£l„/i| eV, where is the Hubble parameter 
in units of 50 km s _1 Mpc _1 , using the current CMB temperature of 2.726 K, and relevant parameters from 
Kolb & Turner (1990). The present group analysis' favored £1„ ~ 0.20 then corresponds to m„ ~ 4.6 eV, for 
H = 50 km s _1 Mpc _1 . 

The void probability function is perhaps another indication that CHDM works better for lower £1„ 
(Ghigna et al. 1994). Also, CHDM forms structure so late that observations of damped Lyman alpha 
systems in quasars are now only marginally consistent with the CHDM parameters studied here (Klypin, 
et al. 1994a and references therein). However, the abundance of massive objects at high redshift is quite 
sensitive to £1„, since these objects are far out on the exponential tail of the mass distrubution. Reducing 
£1„ to 0.20 or even 0.25 brings good agreement with current observations (Klypin, et al. 1994a). 

The second major conclusion of the present paper is that, within the models studied here, it is quite 
natural to find £1 ~ 0.1 from group mass to light ratios, even though the true SI is 1. Three factors combine 
in the same direction to severely bias the M/L method on the low side. First, galaxies within a group 
occupy only the central core of much larger DM concentrations. The total mass bound to the group at 
2r/j is a factor of 3 higher than the DM within r^. Second, the appropriate virial mass to calculate is that 
due to the individual cold DM particles, not the galaxies, whose velocity bias contributes another factor of 
~ 2. Finally, the virial theorem implicitly assumes spherical symmetry, and our groups are actually fairly 
elongated. Counting mass within these elongated boundaries will add perhaps another 30% to the true mass, 
giving a net correction to M/L of nearly an order of magnitude. Hot DM which is unclustered will further 
bias M/L to the low side. 
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Figure Captions 



Figure 1. Comparison of the two methods of assigning velocities to the fragments of (a) CHDM and (b) 
CDM1 overmerger breakups. V°™ gh , the rms velocity of all neighboring galaxies within 1 Mpc , or within 
whatever larger radius encloses 4 neighbors, has a much wider range than the characteristic circular velocity 
V c om = (GM eJJ /r eJj y/ 2 , yet the medians are very similar. CDM1 has a wider range for both velocities than 
does CHDM. 

Figure 2. Comparison of the halo 1-cell mass distribution for the different massive halo overmerger breakup 
methods. Our favored method is labeled "breakup", while Methods 1 and 2 are labelled "breakupi" and 
"breakup2". ALogN / AlogM = — 1.3 below the breakup mass Mi, u . 

Figure 3. The DM density distribution for (a) the CHDM simulation box groups made from halos with 
6p/p > 80, and (b) the CDM1 and CDM1.5 box groups made from halos with 6p/p > 156, compared to that 
for r -2 isothermal spheres. r=0 is the uweighted average of the position vectors of all member halos. Only 
the curve for total mass is shown for CHDMi, while CHDM2 shows the hot and cold distributions as well. 
Even at r = 2r^, mass is substantially above the critical density p c . All simulations show an exponential 
p(r), but CDM shows a steeper fall-off than CHDM. Data inside ~ I9bkpc is unresolved and not plotted. 

Figure 4. The cumulative mass distribution around (a) the fiducial box groups and (b) around the fiducial 
sky groups. There is substantial mass outside but bound to the box groups; M(2 r h)/M(rh) is 2.7 for CHDM 
and nearly 2.0 for CDM. Around the sky groups, M(2 r h)/M(rh) is as high as 3.0 for all simulations. Blurring 
in redshift space reduces the differences between CHDM and CDM curves in the sky groups. Data inside 
~ I9bkpc is unresolved and not plotted. 

Figure 5. The ratio of cold to hot dark matter for the CHDM simulation box groups. This is an unweighted 
average over all groups. It is more representative of small groups, which make up most of this volume- 
limited sample. Group cores are relatively dense and cold, with the hot particles predominately forming a 
more diffuse background. 

Figure 6. (a) The assigned 1-cell mass-to-light ratio required in order to reproduce the CfAl Schechter 
luminosity function parameters, for each simulation. There is only a single curve for the Cen-Ostriker and 
Tully-Fisher prescriptions, since no attempt is made to force the resulting luminosity function to agree with 
observations. The steep rise in M/L at low mass means low mass galaxies do not appear in the sky catalogs, 
(b) The distribution of no breakup sky catalog galaxies vs. mass, (c) Breaking up overmergers clips the high 
end and makes a large number of fragments below Mi, u . Tully-Fisher (not shown) and Cen-Ostriker L(M)'s 
put proportionately much higher luminosity on faint, low mass, nearby galaxies. 

Figure 7. The median axial ratio (depth divided by average of RA and Dec dimensions) in real space of 
redshift-selected groups, vs. V5. The CHDM curves flatten at low V5, indicating minimal contamination 
by interlopers here. CDM curves are steeper and more contaminated. For comparison, Poisson distributed 
galaxies of rms peculiar velocities a = 200 and 350 km s _1 are much steeper still, showing even CDM galaxies 
have significant coherent motion. 
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Figure 8. The number-weighted virial mass estimator M„; r from the halos, vs. the dark matter Mum 
contained within the mean harmonic radius defined by the halos, for 3D selected box groups in CDM1. 
The no-breakup case (a) has a significant tail of high M„; r / 'Mum groups, while the breakup case (b) includes 
many more M„; r / 'Mum ~ 1 groups. Higher pairwise velocities in CDM cause the extended tail to high 
M v i r I 'Mum in both cases. The combined sky groups for all six viewpoints (c) show a strong bias towards 
high M v i r I 'Mum at low mass, as high pairwise velocities mean small sky groups are not bound. 

Figure 9. The same as Figure 9, but for CHDM2. There is little difference between the (a) no-breakup and 
(b) breakup cases, and the low pairwise velocities make for less bias and less noise in the (c) sky groups' 
M v i r than seen in CDM. 

Figure 10. The velocity bias parameter b v for all box groups at the fiducial link. The horizontal lines define 
the median values. CHDM2 has a significant bias of b v = 0.7. CDM1 shows almost no net velocity bias. 
CHDMVs low pairwise velocities lead to a tighter distribution. The banding in CDM1.5 (d) is due to the 
small subsample of dark particles used in the calculations. 

Figure 11. The galaxy density vs. redshift distribution in the magnitude limited sky catalogs for CHDM2 
vs. that for CfAl. Since the simulations' densities drops less steeply with redshift due to the limited box size 
and cosmic variance, requiring random sub-samples of simulation groups to match CfAl's group distribution 
N(z) insures that median properties have no distance-dependent bias. This is the N(z) correction. 

Figure 12. Median v gT vs. V5, for sky catalogs from the breakup boxes and the CfAl. CHDM curves match 
the observations well, while CDM's are too high. Below ~ 400 km s _1 , curves are degenerate since the CDM 
groups are in the "clipped" regime. On this and later curves, error bars are la over the 6 viewpoints. 

Figure 13. f gT vs. V5 for all simulations. CDM groups too few galaxies. CHDMi's higher power on large 
scales leads to grouping too high a fraction, while CHDM2 agrees well with observations. 

Figure 14. Our favored statistic v gT vs. f gT under varying V5 for our fully corrected case using (a) the 
10.384 sr sky catalogs and (b) the 2.66 sr CfA-sky catalogs. CHDM2 is in close agreement with the CfAl 
data. CDM curves in (a) are too high by an estimated ~ 6<r due to cosmic variance, assuming the difference 
between CHDMi and CHDM 2 is ~ 1.7<r. 

Figure 15. Comparison of v gT vs. f gT for the fiducial (a) Schechter, (b) Cen-Ostriker, and (c) Tully Fisher 
luminosity assignment methods. Comparisons were done before breaking up overmergers or applying an N(z) 
correction, and the purpose here is mainly to demonstrate robustness. In all cases, CDM is substantially 
too high. The Tully-Fisher case is probably the least reliable method, as it produces a luminosity function 
which is quite un-Schechter-like. All methods give similar conclusions, though TF predicts a lower £1„ than 
the others. 

Figure 16. Comparison of v gT vs. f gT for different overmerger breakup methods for (a) CHDM2 and (b) 
CDM1. Other simulations are similar. All methods except Method 1 all lead to CDM curves several a 
too high. Method 1 collapses all simulations onto the same curve (e.g. (a) and (b) dotted curves virtually 
overlap) by forcing all fragments to have the maximum possible nearly equal brightness. 
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Figure 17. Median v gT vs. / gr when varying D n and holding V5 constant at the two different fiducial values 
(a) V5 = 600 and (b) V5 = 350 km s _1 . These comparisons were done with no breakup of overmergers and 
no N(z) correction. v gT curves are fairly flat, since groups grow mostly on the sky direction when varying 
D n . At V5 = 350 the CDM groups are in the "clipped" regime. A better comparison is at V5 = 600 km s _1 . 
NKP94 used this to (prematurely) claim CfAl groups had significantly cooler cores than any simulations. 

Figure 18. The same as Figure 17, but using the breakup catalogs and making the N(z) correction to the 
median v gT 's. The "cooler cores" conclusion of NKP94 is now seen to be an artifact of overmerging and 
N(z) bias. When corrected, CHDM is in good agreement with observations at both fiducial Vs's, while CDM 
remains too high. 

Figure 19. The same as Figure 18, but now using the CfA-sky catalogs. The error bars are larger, but 
otherwise the figure is virtually identical to that for the full 10.384 sr catalogs. At low and moderate V5, 
/ gr differs very little between these samples. 

Figure 20. Percolation properties of groups vs. varying V5 in the fully corrected case for the CfA-sky catalogs. 
CDM produces too few groups at low V5, while CHDM2 produces too many at intermediate and high V5. 
CHDMi percolates too easily, as groups merge faster than they are created above V5 = 450 km s _1 . A curve 
intermediate between CHDMi and CHDM2 would likely give the best fit, though discrimination between 
models is poor in this plane. 

Figure 21. Percolation properties of groups vs. varying D n for all simulations, using the CfA-sky catalogs. 
As in Figure 20, CHDMi percolates too easily. CHDM2 percolates correctly, but produces too many groups 
at all D n . A CHDM model with large scale power intermediate between these, and lower £1„ would appear 
to fit better. Both CDM models fit well on this measure. 

Figure 22. The ratio of the simulation mean harmonic radius to that for the CfAl at the same links, 
vs. D n . All simulations are significantly too high. This may be due to residual resolution limitations, or to 
simulation halos having spatial distributions not sufficiently galaxy-like. 
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